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We propose a novel interaction-based route to half-metal state for interacting electrons on two-dimensional 
lattices. Magnetic field applied parallel to the lattice is used to tune one of the spin densities to a particular 
commensurate with the lattice value in which the system spontaneously 'locks in' via van Hove enhanced 
density wave state. Electrons of opposite spin polarization retain their metallic character and provide for the 
half-metal state which, in addition, supports magnetization plateau in a finite interval of external magnetic field. 
Similar half-metal state is realized in the finite-U Hubbard model on a triangular lattice at 1/3 of the maximum 
magnetization. 
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Spin systems supporting robust magnetization plateaux 
whereby macroscopic magnetization M is fixed at a rational 
fraction of the full (saturated) magnetization value M sat in a 
finite interval of external magnetic field hi < h < hi are 
subject of intense experimental studies [1-5]. Typically these 
materials are magnetic insulators which are well described by 
the Heisenberg-type models with short-range exchange inter- 
actions between localized spins. 

One of the best understood and studied plateau states is rep- 
resented by the up-up-down (UUD) magnetization plateau at 
M = j M sat in the triangular lattice antiferromagnet [6, 7]. 
This 1/3 magnetization plateau is a remarkably stable state - 
it is known to survive significant spatial deformation of ex- 
change integrals in both quantum (spin 1/2) and classical ver- 
sions of the model [8, 9], well beyond the point where adjacent 
to it co-planar spin states ceases to exist altogether. The basic 
reason for this stability lies in the collinear structure of the 
UUD configuration. Collinearity preserves U(l) symmetry of 
the Hamiltonian with respect to the magnetic field axis. The 
only symmetry that the UUD state breaks is the discrete trans- 
lational symmetry - its unit cell consists of two up and one 
down spin. This ensures the absence of the gapless (Gold- 
stone) modes in the spectrum and implies enhanced stability 
of the collinear spin arrangement. 

Since the Heisenberg model is just a low-energy approxi- 
mation to the large-[//t limit (here t is the hopping integral 
and U is on-site interaction energy, see below) of the half- 
filled Hubbard model, the insulating magnetization plateau 
state is favored by strong electron-electron interactions. What 
happens to the 1 /3 magnetization plateau state as the ratio U/t 
is reduced and electrons delocalize is one of the key questions 
motivating our study. 

A different class of magnetization plateau materials is pro- 
vided by half-metallic ferromagnets in which by virtue of pe- 
culiar electronic structure all conduction electrons have the 
same spin orientation (say, up, for definiteness). In their sim- 
plest version half-metallic materials are then fully saturated, 
M = M sa t. As the name suggests, these materials are con- 



ductors and are well understood in terms of non-interacting 
electron picture [10, 11]. 

The aim of our work is to unite these phenomena by propos- 
ing two new interacting routes to the half-metallic magnetiza- 
tion plateau states. Both routes require finite external (Zee- 
man) magnetic field, applied parallel to the two-dimensional 
triangular lattice. 

The weak-coupling route, described first below, relies on 
tuning density of majority (say, spin-up) electrons n-j- to a spe- 
cific value (3/4), commensurate with the triangular lattice, at 
which the Fermi surface (FS) passes via a set of van Hove 
points with logarithmically divergent density of states, see 
Figure 1. Depending on the total electron density n = rif+n^, 
the FS of minority (spin-down) electrons may or may not be 
affected by the interactions, but in any case retains its metal- 
lic character. The resulting ground state is half-metal which 
supports M = (| — n^)M sa t magnetization plateau with ferri- 
magnetic (up-down-down-down) collinear spin structure. We 
emphasize that M/M s t a t is generally irrational. This novel 
state has no analogs in the in the large- U limit of the Hubbard 
model and can be pictured as a phase with coexisting spin- 
and c/iarge-density wave orders. Theoretical analysis of this 
limit bears strong similarities with recent proposals [12-15] 
of collinear and chiral spin-density wave (SDW) states of itin- 
erant electrons on honeycomb lattice in vicinity of electron 
filling factors n = 3/8 and 5/8 and at zero magnetization. In 
contrast to our problem, however, SDW order there sponta- 
neously breaks spin-rotational symmetry of the Hamiltonian 
and, as a result, is accompanied by gapless collective excita- 
tions [16] which drive the competition between the collinear 
and chiral orders at finite temperature [13]. This complication 
is absent in our problem where external magnetic field sets di- 
rection of the collinear SDW. The resulting half-metallic state 
breaks only discrete translational symmetry of the lattice and 
is stable with respect to fluctuations of the order parameter 
about its mean-field value. 

Next we describe the strong coupling (large-C/) route 
and show that M = jA/ sa t magnetization plateau (n-f = 



2 



ky 



- 










: / 


(J— 


V 02, 


QiV 














j-j 

\Q3 
















FIG. 1 : (Color online) The non-interacting Fermi surfaces of spin-up 
and spin-down electrons (top) at 1/2 magnetization and the recon- 
structed fermi surface of spin-down electrons at Vi/t = 0.1 (bot- 
tom). The thick black hexagon is the first Brillouin zone of triangular 
lattice. The red hexagon is the Fermi surface of the spin-up electron. 
It is perfectly nested by linear combinations of three wave vectors 
(the arrows). The dashed hexagon is the first Brillouin zone under 
the folded scheme. The dashed purple circle is an example of the 
hot-spot-free Fermi surface of spin-down electrons in the hole doped 
system. 



2/3, ni = 1/3), present in the U/t — > oo limit, survives 
down to Ud/t w 4.3, which is significantly lower than the 
zero-magnetization critical value t/120 /t ~ 10 below which 
the three-sublattice 120° magnetic order melts as the system 
transitions to a quantum spin-liquid state [17-19]. Quite in- 
terestingly, we find that in the interval U c \ < U < U C 2 ~ 4.8i 
the UUD state is a half-metal with mobile majority (spin-up) 
electrons. A closely related question of the evolution of the 
ground state of the Hubbard model at zero magnetization, as a 
function of U/t, has been actively investigated [17-21], partly 
because of the potential relevance to intriguing 'spin-liquid' 
organic insulators [22-24]. 

Weak-coupling analysis is based on the extended Hubbard 



model description of electrons on triangular lattice: 

H = —t (c\ a C r t a + ll.c) + U 22 n rt n rl + V n r n r' 

(rr 1 ) r (rr') 
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Here (rr') labels nearest neighbour bonds, < U/t, V/t <C 1 
are onsite and nearest neighbour repulsive interactions, re- 
spectively. For now we set the filling factor at n = (n/) = 
(ntf + riii) = 1. 

The Zeeman field h, normalized to include the usual qhb 
factor, is tuned to produce M = ^M sat magnetization so that 
on average n-j- = 3/4 and = 1/4 per site. Under this condi- 
tion the Fermi surface of spin-up electrons, shown in Figure 1, 
is given by a perfect hexagon whose vertices, located at the M 
points of the first Brillouin zone (BZ), are van Hove (saddle) 
points of the dispersion with vanishing Fermi velocity [13]. 
These saddle points are the reason for the (logarithmically) di- 
verging density of states as well as singular susceptibility (see 
below). They are connected by the wave vectors Qi = ^=y 
and Q2,3 = T^x — -^y (Fig. 1), which are just halfs of the 
corresponding reciprocal lattice vectors G 1,2.3 [12]. In ad- 
dition, parallel faces of the FS are perfectly nested by linear 
combinations of Q's. Spin-down FS is nearly circular, as Fig. 
1 shows, and does not possess any special features. 

The very special role of the wave vectors Qi,2,3 is con- 
veniently quantified by charge susceptibility X CT (q) of spin-er 
electrons defined in the standard way as 
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Here N is the number of sites, rik., is the occupation number 
of fermions with spin a and momentum k and 
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is the free-particle dispersion. Straightforward calculation, 
outlined in [25], shows that susceptibility of spin-up electrons 
is strongly divergent at q = ±Q a , while that of spin-down 
electrons is finite, 

*t(Q.) - -^ 2 (^t = ^, (4) 

Xl(Q a ) = C i w 0.136. (5) 

Here A ~ ir/a is the 'size' of the BZ (a is the lattice spacing) 
while q ~ 1/L is the microscopic cut-off which scales as the 
inverse linear size of the lattice L ~ V^A a. 

Eq. (4) suggests strong modulation of density at q = ±Q a 
which motivates the following mean-field ansatz 
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where index a describes two spin projections, a =t= +1 and 
a =1= — 1, in the usual way. Note that the filling factors 
n a = J2 r (n r a ) /N are not affected by finite order parameters 
m a , and = 3/4, while = 1/4. Ansatz (6) allows us to 
approximate ( 1) as 



k,cj k,(T,a 

+3NV(m t + m l ) 2 ~ ZNUm^m^ 
y a = Um_ a - 2V(m t + mj.), 



(7) 
(8) 



The amplitudes of density modulation of spin-cr electrons are 
determined self-consistently 



— ^cos(Q a • r)(n r>0 
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where dispersion is determined by the mean-field Hamil- 
tonian (7). We find = Ci(V — V /2)m\/t, which means 
1 7714.1 < |m T |, and 



(10) 



The main effect of interaction y-j- in (7) is to provide an infra- 
red cut-off q-f ~ l^/il 1 / 2 in the susceptibility (4) of spin-up 
electrons at the van Hove points. This leads to the final result 
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(11) 



Non-analytic dependence of m on interaction amplitudes U, V 
is determined by van Hove points. The sign of m-j- is chosen 
such that the FS of spin-up electrons is gapped for all mo- 
menta. We checked that the opposite sign leads to a quadratic 
touching of the top two bands at T point and results in a state 
of higher energy [25]. 

Equations (6) and (11) show that the ground state is a super- 
position of commensurate charge-density and collinear spin- 
density waves. This is a four-sublattice state: observe that 
J2 a =i C0S (Qa • r ) takes values 3, —1, —1, —1 on the sites of 
triangular lattice r = (x,y) = diai+<i2a2, where ai = (1,0) 
and a.2 = (1/2, \/3/2) are elementary lattice vectors and 
are integers. We stress that the perfectly nested FS of spin-up 
electrons is crucial for the spin-up electrons to be gapped at 
arbitrary weak interaction. 

Note that in the absence of direct density-density inter- 
action, when V = 0, the order parameter scales as m ~ 
cxp (—const/ U). The density wave of spin-up electrons in 
this case is driven by the effective interaction cx X1U 2 , which 
is mediated by spin-down electrons, since the onsite repulsion 



U only couples electrons of opposite spins. A small finite 
V > U 2 /t, see (11), changes this scaling to a much stronger 

dependence m ~ exp (— const/ \fv\ . 

The band structure of spin-down electrons is also modified 
as (7) shows. For finite 7/4, while the spin-down electrons 
remain gapless, the 'hot spots' on the spin-down FS (which 
are the points connected by Q Q ) are gapped and the FS is re- 
constructed as shown in Figure 1. By reducing the density 
of spin-down electrons below 1/4, while maintaining that 
of spin-up ones at the perfect nesting condition = 3/4, 
one can reduce the spin-down FS below the critical volume to 
make it fit inside the reduced Brillouin zone (dashed hexagon 
in Fig. 1), which is 4 times smaller than the original one, 
and the hot spots disappear altogether. We find that this hap- 
pens for n cr < 0.976 (Figure 1), i.e. n CT i < 0.226. Under 
this condition, and for weak interactions < t, the FS 
of spin-down electrons is not affected by the reconstruction of 
the FS of spin-up electrons at all. In either case, the result is a 
half-metal where all conducting electrons have spin opposite 
to the direction of the external field. 

Similar half-metallic state can be realized on other lattices. 
For example, consider the same Hamiltonian (1) on the square 
lattice with filling factor n ^ 1. A proper field h can be ap- 
plied such that the FS of the spin-up (if total density n < 1) or 
the spin-down (if n > 1) electrons is perfectly nested. Clearly 
this FS is a square with vertices at (±7r, 0) and (0, ±7r) points. 
The system then develops density wave order at wave vector 
Q = (71", 71"). The ordered state is generically half-metallic. 
Note also that in this case the FS of conduction electrons ex- 
periences no reconstruction due to the absence of hot spots. 

Strong coupling limit. We now turn to the question of 
M = |M sa t plateau which exists in the opposite limit of 
strong interactions, U ^> t. As should be clear from the previ- 
ous discussion, in the current situation filling factors of spin- 
up (down) electrons 2/3(1/3) do not correspond to perfectly 
nested FS. In order to make connection with the insulating 
magnetization plateau phase of the Heisenberg spin model we 
set V = 0, keep total density at one electron per site n = 1 
and introduce the following mean-field ansatz: 
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— 2rj a cos(Q • r) 



(12) 



where Q = ^-x describes the UUD pattern. (Note that 
cos(Q • r) takes values 1,-1/2,-1/2 on the triangular 
lattice.) Parameters i] a of the state are determined self- 
consistently by the equations similar to (9). Solving them nu- 
merically we find discontinuous jump of r\ a from zero to finite 
values when U > 4.30i for a range of h. We also find that rj-f 
and ifi are l n general different so that the system displays a 
co-existence of the spin-density and charge-density wave or- 
ders [25]. Interestingly, it is now the spin-down electrons that 
are gapped while the spin-up electrons remain gapless around 
the Fermi energy, see Figure 2. Spin-down electrons fill com- 
pletely the lowest of the three bands wj,(k) in the folded Bril- 
louin zone, while the spin-up ones fill the two lowest bands 
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FIG. 2: (Color online) The spin-up (top) and spin-down (bottom) 
bands of the half -metal state at 1/3 magnetization at U = 4.60t. 
In order to give an idea of dispersion along the k y axis as well, 
the graph shows projections of two-dimensional bands onto k x axis 
for 30 different discrete k y values. The Fermi energies fi a — 
(i + a(h/2 + U /6) are shown as black-dash lines. 



w^(k). For even stronger interaction U > 4.8CM, the two up- 
per spin-up bands also get separated by a gap which turns the 
half-metal state into an insulator with collinear UUD pattern 
of local magnetization. 

To estimate how competitive the found half-metal state is 
we compared its energy with the uniformly magnetized trans- 
verse spin-density wave state, i.e. "cone" state in magnetic 
language. This state is characterized by the two order param- 
eters, longitudinal (S% ) = M and transverse = moe ,q ' r 
magnetizations [26, 27], where the ordering wave vector q is 
incommensurate in general. Comparing their mean-field en- 
ergies, we determined that the half-metal state has a lower 
energy for 4.45£ < U < 4.80t [25]. The resulting mean- 
field phase diagram is shown in Figure 3. Note that both half- 
metal and insulator phases are magnetization plateau states 
with M = ±M sat . 
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FIG. 3: (Color online) The mean-field phase boundaries of the half- 
metal and insulating states at M = 1/3 M sat on triangular lattice. 
Solid (dashed) lines denote first (second) order transitions. 



Discussion: we described two general ways to induce half- 
metal states through a combination of electronic interactions 
and a finite Zeeman field. Our work provides new avenue 
to half-metallic states which have potential applications in 
future spin-dependent electronics. Despite our focus on the 
Hubbard-Heisenberg model on the triangular lattice, both of 
the proposed mechanisms can be easily generalized to other 
lattice geometries. It is important to note that in our proposal 
the Zeeman field sets direction of the spin polarization of con- 
ducting electrons in the half-metal state. This convenient fea- 
ture makes it different from the recently discussed half-metal 
SDW state at zero magnetization [13] where spin-rotational 
symmetry breaks spontaneously, resulting in a random selec- 
tion of the spin polarization axis in the spin space. 

Our proposal can also be realized in Kondo lattice systems 
where the Zeeman field is provided by exchange couplings 
between the itinerant electrons and ferromagnetically ordered 
local moments, see e.g. [12], as well as in cold atom systems 
[28, 29] where it may be easier to achieve the required spin 
population imbalance. 

Magnetic field control of the half-metallic phase may be 
useful for creating switchable interfaces between half-metal 
and noncentrosymmetric superconductor which have been ar- 
gued to support Majorana bound states [30]. 

It is worth noting close physical similarity between our pro- 
posal and the previously proposed one-dimensional 'Coulomb 
drag' setup [31] where role of the lattice is played by the elec- 
trons in an active wire interactions with which gap out one of 
the spin projections in the passive wire. 

Several interesting theoretical questions can be asked re- 
garding the half-metal state. First of all, the metal-insulator 
transition between a half-metal and a Mott insulator, found 
here at U/t w 4.8, represents Mott transition which is not 
affected by (gapped) spin fluctuations. Understanding it in 
details may lead one to a better characterization of the general 
Mott transition. Half metal states are adjacent to many other 
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interesting quantum phases. For example, it was demonstrated 
[32] that d+id chiral superconducting state is the ground state 
at 3/2 (5/8) filling on triangular (honeycomb) lattice for weak 
interaction. If we dope the system with holes while keeping 
the FS of spin-up electrons perfectly nested by a Zeeman field, 
the system will eventually become a half-metal. Understand- 
ing how quantum phase transition(s) between these different 
phases happen is an interesting question left for future studies. 
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like to thank for hospitality. The work is supported by NSF 
through Grant No. DMR-0808842 (O.A.S.) and by NSERC 
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Supplementary Material for "Half-metallic magnetization 
plateaux" by Zhihao Hao and Oleg A. Starykh 



Susceptibility 

Susceptibility Xo-(Qo) (a = 1,2,3) of free electrons is 
given by 



XAQa) 



2, n k,cr — «k+Q a ,<r 



d z k 



(1) 



where the integration is restricted in the first Brillouin zone of 



2 cos cos 



area A = 8ir 2 /y/3. = — 2 ^cosfc^ 
is the dispersion of the non-interacting fermion. Due to six- 
fold rotational symmetry, x CT (Qa) is independent of a. We 
consider a =f and a = 1. The particle (n^ = 1) and hole 
(?ik = 0) states contribute equally to x- in addition, the 
integrand is invariant under separate operations k x — > —k x 
and k y — > —k y . Taking advantage of these symmetries, we 
compute the leading singular contribution to x by writing 
k = (tt, + q and considering small q limit: 



Xt(Qi 
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dq x 
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For a =4-, Xi is computed numerically: 

0.136 
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(3) 



Half -metal state at weak coupling 




FIG. 1: The Fermi surface at = (the blue lines) under the 
folded zone scheme. The folded zone (the red hexagons) is four 
times smaller than the first Brillouin zone of the triangular lattice 
(the black hexagon). 
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We consider only spin-up fermions on triangular lattice at 
3/4 filling. The following Hamiltonian is studied: 

H = H U ek ~ ^) c kt c kt + | X^( c kt Ck +Q 5 t + h - c ) • 

k V a=l / 

(4) 

The Hamiltonian can be diagonalized and we obtain four 
bands in the folded zone scheme. In the limit of y —> 0, three 
of the four bands become £k+Q a (a = 1, 2, 3). These bands 
touch at three straight lines through the V point (Fig. 1). The 
spectrum at k = is —8, —y, —y and 2y. For y > 0, the 
lower three bands are —8, — y and —y. The Fermi surface is 
fully gapped. For y < 0, the Fermi surface become a point at 
k = 0. 

We define the energy of the Fermi sea to be E(y). At small 
\y\, the band structures of y > and y < are only signifi- 
cantly different around the T point. E(y) for positive y is low- 
ered by order \y\ in an area proportional to |y| around k = 
since the energies at the van Hove point are functions of k 2 for 
small k. SE(\y\) = E(—\y\) — E(\y\) can then be estimated 



FIG. 2: 5E(y) / In 2 (j/ 1 ^ 2 ) (blue squares) as a function of y fitted by 
function 0.055y 2 (the red line). 



as proportional to Xt(M)y 2 ~ l n (\y\ 1 ^ 2 )y 2 > as shown by 
the red line in Fig. 2. The blue points are obtained by numer- 
ically calculating 5E(\y\) as a function of momentum k and 
summing over fine mesh of fc-points in the Brillouin zone. 



M = M sat /3 half -metal state 
We use the following mean-field ansatz: 

(n rCT ) = - + ^ - 2rj a cos(Q • r). 
2 o 

The mean-field Hamiltonian can be written as: 



2 




form magnetization. We use the following notation for order 



FIG. 3: The Fermi surfaces of spin-up (blue) and spin-down (purple) 
electrons at M = | M sa t for r) a = 0. 
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FIG. 5: Qi mm V(47r/3) as a function of U 




FIG. 4: r\ a as a function of U. 
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The self-consistent equations read: 

^ = _ ^5Z C0S ^ a ' r )(4a c r<T>, (6a) 



1 = ^E c ™<w- (6b) 



These equations are solved numerically. r\ a become finite for 
U > AM (Fig. 4). 

We compare the energy of the half-metal state with a sec- 
ond mean-field ansatz: a spiral spin-density wave with a uni- 



Define the two component spinor as V'k = (ck+Qf 5 Ck^) T , 
the Hamiltonian can be written as: 

H = ^^(e k+Q + e k - / x)l4-(e k+Q -e k -~-tfM)a 2 
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V'k + NUM 2 + NUml. 
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where e k = — 2t(cos(k x ) + 2 cos(k x /2) cos(VSk y / 2)) . The 
spectrum can be solved: 
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The self consistent equations read: 
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Here 0(x) = 1 for x > and 8 (a:) = for x < 0. 

For fixed U and /i, we solve the self-consistency equations 
numerically for different Q's. Focusing on Q = (Q x ,0), we 
find the Q( min ) with lowest mean-field energy. Q( mm ) is in- 
commensurate in general (Fig. 5). We scan over the two di- 
mensional phase space Q for some U's and confirm that the 
lowest mean-field energy state is indeed on the Q y = axis. 

Comparing the mean-field energy of the two states, the half- 
metal state has lower energy for U > 4.45t for ranges of h. 
For U > 4.81 the half-metal state becomes a band insulator. 



